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! Abstract 



We study nonleptonic A& — ► A c 7r, S c 7r and Y,*ir decays in the limit m&, m c , E n 3> Aqcd using 
the soft-collinear effective theory. Here S c = S c (2455) and S* = S c (2520). At leading order the 
^p" 1 . Af, — ► si*^7r rates vanish, while the A& — ► A c 7r rate is related to A& — > A c lv, and is expected to be 

larger than r(£> — > D^tt). The dominant contributions to the A& — ► e£^7t rates are suppressed by 



A qcd/^- We predict T(A 6 -> S*vr)/r(A b -► S c vr) = 2 + [Aqcd A^Q , a,(mg)], and the same 



^ ■ ratio for A b -» E£> and for A b -> H^if. "Bow tie" dia grams are shown to be suppressed. We 

comment on possible discovery channels for weakly decaying pentaquarks, 0{, iC and their nearby 
heavy quark spin symmetry partners, 0£ c . 
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Heavy baryon decays are interesting for many reasons. Heavy quark symmetry ^] is more 
predictive in semileptonic A b — > A c tv decay than in B — > D^Hu, and may eventually give a 
determination of jV^il competitive with meson decays Q|. In this paper we concentrate on the 
more complicated case of nonleptonic b —>■ cud baryon transitions, as shown in TableUJ These 
channels provide a testing ground for our understanding of QCD in nonleptonic decays. Our 
analysis is based on heavy quark symmetry and the soft-collinear effective theory (SCET) 

There is considerable experimental interest in these decays. Recently the CDF Collabo- 
ration measured 

f Ah B(A b ^A+7T-) 

U W° - D + n-) = °- 66 ± °- 11(stat) ± °-° 9(syst) ± °- 18(BR) ' (1) 
where f\ b and fd are the fragmentation fractions of b quarks to A b and B°, respectively. 
Using the input f bary(m / f d = 0.304 ± 0.053, CDF quoted B(A b -> A+ir-)/B(B° -> D + ^-) ~ 
2.2 The measured lifetimes, r(A&) = 1.23 ps and t(B°) = 1.54 ps p, then imply that 
r(Afe — > A^tt~)/T(B° — > D + 7i~) ~ 2.7. More experimental results on semileptonic and 
other nonleptonic A fe decays are expected in the near future. 
The part of the weak Hamiltonian relevant for this paper is 

4G F 



H w - VcbV* d 



C 1 (m b )O l (m b )+C 2 {m b )0 2 (m b ) , (2) 
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TABLE I: The decays considered in this paper. Here si is the spin of the light degrees of freedom 
0. The mass shown for the S c *^ is the average of the charge and +1 states, and for the S c 's the 
mass is the average in the doublet. The stability of the pentaquark states 0q(= Qudud) and their 
value of S[ are both conjectures. 
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Tree Color-commensurate Exchange Bow tie 

FIG. 1: Diagrams for A& decays, giving amplitudes T, C, E, and B. Decay to A c gets contributions 
from all four terms. Decays to Tr c and H c do not have T and T, C contributions, respectively. 



where both the Wilson coefficients, Cj, and the four-quark operators 

Ox = (c^P L b) (d>fP L u) , 2 = (d^P L b) {cYPlu) , (3) 

depend on the renormalization scale which we take to be m&, and Pl = (1 — 7s)/2. The 
combination [Ci(m&) + C2(m&)/3] is very close to unity. 

Weak nonleptonic decays are sometimes characterized by diagrams corresponding to 
different Wick contractions. As shown in Fig. ^ there are more possibilities in baryon 
than in meson decays. In particular, a "Bow tie" contraction is unique to baryons. The 
color structure for baryons also differs from mesons: we find that the C diagram is of 
the same order in the large N c limit as the T diagram. 1 Nonleptonic meson decay ampli- 
tudes are sometimes estimated using naive factorization, which would set (A c 7r|Oi|A&) = 
(A c |c7 m Pl6|A{,) (ir\d'~f fl PLu\0) . In baryon decays the extra light quark implies that this pro- 
cedure is ill-defined for all but the tree diagram. In naive factorization the A;, — > T,^ir 
decays are very suppressed, since the T contribution vanishes separately in the isospin and 
heavy quark limits |7[ (just like the semileptonic A& — > E[*'^ decays), the C contribution 
vanishes after doing a Fiertz transformation on the four-quark operator, and the E and B 
amplitudes are identically zero since the u and b fields are in different quark bilinears. 

In this letter we show that more rigorous techniques can still be applied to make rea- 
sonable predictions for all these decays. By expanding in m&, m c , E n ^> Aqcd we show 
that for Ak — > A^tt" the amplitudes corresponding to the diagrams in Fig. ^ satisfy 
T ^> C ~ E ^> B, and we find that the experimental result in Eq. (JTJ) is consistent with 

1 If we treated the N c — 3 additional quarks in the baryons as flavors that arc sterile under the weak 
interaction then color-commensurate would become color-suppressed. 
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theoretical expectations. Next we consider A b — > S^7r decays, and show the leading con- 
tributions to these nonleptonic rates are suppressed by Aq CD /E%, much like in B° — > D°n°. 
Using heavy quark symmetry we derive a relation between the decay rate to S c and £* and 
comment on decays to S c . Finally we consider the detection of possible weakly decaying 
heavy pentaquarks, Q b and C , with nonleptonic decays. 

The proof of factorization at leading order for A b — > A c tt decay follows closely that for 
so we do not review it here. In this case the nonperturbative expansion 
parameter for SCET is A = Aqq^/E^ Since E n is set by the bottom and charm quark 
masses, we take this to be of the same order as the expansion parameter for the heavy quark 
effective theory (HQET), i.e., A ~ AqcdA^q (Q — c,b). Working at leading order in A 
and a s (m b ) and neglecting the pion mass, the A b — > A c tt matrix element factorizes in the 
standard way, 

(A c (v', s') n\H w \A b (v, s)) = v^C?f(Ci + ^V cb V: d f n E n (A c (v f , s')\c^P L b \A b (v, s)) , (4) 

where f n = 131 MeV is the pion decay constant, n is a light-like four-vector along the 
direction of the pion's four-momentum, p% = E^^, and the four- velocities of the A b and A c 
are v and v', respectively. Perturbative corrections induce a multiplicative factor in Eq. 
(T(x)) 7r = Jq dxT(x) 07r(x), where T(x) is computable and (f> n is the nonperturbative light- 
cone pion distribution function and a term proportional to the matrix element of 
cifiPfib. At leading order in a s (mg), we can set (T(x)) n = 1 and the term involving c^P^b 
to 0. This implies that the nonleptonic rate is related to the semileptonic differential decay 
rate at maximal recoil, 

3ir 2 {d + C 2 /3) 2 \V ud \ 2 f 2 ( &T{A b ^A c lv)' 
m\ ta y dw 

where r\ = m\Jm\ b , w = v ■ v ' = (m\ b + m\ c — q 2 )/ (2m^ b m Ac ), and io max corresponds to 

q 2 = m£(~ 0). 

The semileptonic A b — * A c iv form factors are 

(A c (p', s')\V^\A b {p, s)) = u(p', s') [f llfl + f 2 v^ + hv'^uip, s) , 

(A c (p / ,s')|A /Lt |A 6 (p,s)) = U(p', S')[^l7/. + 92Vp + 03^75 «(P, s) , (6) 

where the fa and are functions of w, and the relevant currents are V v = c^i v b and A v = 
crfv^b. The spinors are normalized to u(p, s)7 M w(j», s) = 2p fl . In the heavy quark limit, 

CO) = /iH = 9i(w) , 



= h(w) = f 3 (w) = g 2 (w) = g 3 {w) 



(7) 



where ((w) is the Isgur-Wise function for ground state baryons. The differential decay rate 
is given by 



dT(A b ^AJu) _ G 2 F mj b \V^ 2 

d^ " 24^ ' AV 



3 Vw 2 — 1 Qw + 6wr\ — Ar A — 8r A u> 2 



Jfcw) , (8) 



where in the tuq ^> Aqcd limit T A {w) is equal to the Isgur-Wise function, £(«;), an d i n 
particular ^a(I) = 1. In terms of the original form factors 
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(9) 



F A (w) 2 = [6w + 6wr 2 A - 4r A - 8r A w 2 \ |(iu - 1) [(1 + r A )/i + {w + l)(r A /a + /a) 
+ (w + l) [(1 - r A )#i - (w - l)(r A g 2 + # 3 ) 
+ 2(1 - 2r A u> + r\) \{w - l)f 2 + (w + l)^ 2 



Combining the above results for 
D(* )+ 7r~ we find that 



A+7r decay with the analogous ones for B c 



T(A b -> A c tt-) 



5ml(l-ri) 3 r DW 



A ) 



(10) 



r(jB o _ DW+7r -) ^(1-^)3(1 + ^)2 ve«i; , )y 

where £ is the Isgur-Wise function for B — > semileptonic decay, and r D (*) = m D (*) /mg. 
When the — > A+i'v rate is measured, one can directly test factorization using Eq. (jSJ) or 
Eq. ()10jl. In the absence of this data, we have to resort to using model predictions for the 
baryon Isgur-Wise function. If the ratio of Isgur-Wise functions in Eq. (JTUJl is unity then 
the prefactor in Eq. flTUJ) implies that T(A b -> A c 7r-)/r(5° -> DM+7T-) = 1.6(1.8). This 
enhancement is in rough agreement with the data in Eq. (jlj. A similar result also follows 
from the small velocity limit (rriQ 3> mj — m c ^> Aqcd); i n which the nonleptonic rates 
satisfy T(A b -> A c tt) : r(£ -> £>*7r) : T(B -> Dtt) = 2:1:1, while for the semileptonic 
rates r(A 6 -> A c £p) : r(B -> D*£u) : r(B -> Dfl?) =4:3:1. 

The large N c limit provides some support for the ratio of baryon to meson Isgur-Wise 
functions being close to unity at maximal recoil. In the large N c limit the heavy baryons can 
be treated as bound states of chiral solitons and mesons containing a heavy quark. In this 
picture, the baryon Isgur-Wise function, ((w), is predicted to be ((w) = 0.99 e~ L3< ^ _1 ) jl^. 2 



2 Updating the parameters by fitting the mass splitting to give k = (0.411 GeV) 3 , and using tun = A = 
0.8 GcV (instead of Mjv) for the mass of the light degrees of freedom leaves the exponent essentially 
unchanged. 
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This gives Climax = 1-4) = 0.57, which is indeed close to £(WmIx = 1-5) — 0.55 [3]. Using 
this model for ((w), \V c b\ = 0.04, T\ b = 1.23 ps, and Eqs. © and (JSJ) yield the prediction that 
B(A b -> A+7T") = 4.6 xlO" 3 . As expected, this is larger than 23(1?° -> DW+tt") ~ 2.7xl0~ 3 . 
However, the uncertainty in this prediction is quite large, particularly given that large N c 
strictly only applies for w near 1. The same large N c inputs predict B(A b — > A^£"u) m 6%, 
i.e., this channel is expected to make up a large part of the inclusive A b — > X c £~u rate, with 
the sf = 1~ excited A c states making up a significant fraction of the remainder [l^ . 

Order Aq CD /mQ corrections to these predictions may be significant. The A b —> A+7T - 
amplitude receives contributions from the T, C, E, and B classes of diagrams in Fig. d In 
SCET, \E/T\ and \C/T\ are of order Aqcd/^q |l5j ]. and we will show later that \B/T\ is fur- 
ther suppressed. In B — > Dn decay we know from B(B~ — > D°7T~) / B(B° — > D + ix~) ~ 1.8 0] 
that Aqcd/^q corrections affect the amplitudes at the 15 — 30% level. In particu- 
lar \A(B° -> £) + tt-)| = |T + £| = (5.9 ± 0.3) x 10~ 7 GeV and -> D ^)! = 
\T + C\ = (7.7 ± 0.3) x 10~ 7 GeV. The ratio of these amplitudes can be reproduced 
by a power correction involving a hadronic parameter |s e ff| — 430 MeV, which is of nat- 
ural size Q]. Since B s — > D~ix + only has a T contribution, accurate measurement of 
this rate will improve our understanding of the size of E and C. CDF recently measured 
[f s B(B s - D;ir + )}/lf d B(B° D--K+)} = 0.35 ± 0.05 {stat) ± 0.04 (syst) ± 0.09 (BR ) Q, and 
using f s /f d = 0.26 ±0.03 yields B(B S -> D-n+)/B(B° -> L>-tt+) ~ 1.35. Neglecting 5C/(3) 
breaking 3 this implies \A(B S -> £>Jtt + )| = |T| = (7.3 ± 1.5) x 10~ 7 GeV and that |C| and 
\E\ may be comparable. The errors are still too large to draw any definite conclusions. 

Now we turn to A b — > E c 7r decays. As shown in Table HI there are two E c states with 
different spin which we refer to as E c and E*. They form a heavy quark spin symmetry 
doublet with the spin and parity of the light degrees of freedom, sj' = 1 + . Under isospin, 
the Ah is / = 0, the is J = 1, and the Hamiltonian is I = 1, so the S^7r final state 
must be I = 1 (it can not be I — or 2). Therefore the rates to the two different charge 
channels are equal, 

T(A b sW°7r°) = T(A b £« + . (11) 

Based on 1? decay data and the SCET power counting, we expect T(A b -> E^tt) to be up 
to about an order of magnitude smaller than T(A b — > A c 7r), since the leading contributions 



In the heavy quark limit of the T amplitude ST/ (3) will be tested by the measurement of B s — ► D s tv. 
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FIG. 2: Contributions in SCETi to A b S+vr" [(a) and (b)], and to £°tt [(c) and (d)]. Solid 
dots denote insertions of the suppressed usoft-collinear Lagrangian, C^ q , the double lines are heavy 
quarks, the dashed lines are collinear quarks, the solid lines are usoft quarks, and the "looped lines" 
are collinear gluons. The nonleptonic weak vertex is denoted by ®. 

to Af, — > 'Eyir are power suppressed. 

Again, we use SCET to expand in Aq CD /i71q, Aqcb/E^, and a s (mg), keeping only the 
leading terms that cause the A b — > TS*>it transitions. These come from the C and E diagrams 



in Fig. 
in Ref. 



Q] an 

3- 



and their contributions can be studied following the analysis of B° — > D^n 
The leading diagrams in SCETi that determine the matching onto power 
suppressed operators are shown in Fig. |5J To match the C and E diagrams, two insertions 
of the mixed usoft-collinear Lagrangian, C^) [l^, is required, each yielding a suppression 



of J Aqcb/ E n . This yields the power counting that \C/T\ and \E/T\ are C(Aq CD / E w ). In 
contrast, matching the B diagram in Fig. requires four insertions of (or other higher 
dimensional terms in the Lagrangian), and B is therefore power suppressed compared to C 
and E by at least an additional Aqqd/E^. 

In addition there is a further matching onto SCETn. The resulting matrix element in- 
volves soft and collinear operators which factor Q]. 4 The matrix element of the weak 
Hamiltonian, (£°(i/, s')7r° \Hw\Ab(v, s)), can be written (neglecting Q s (mg) corrections) 
as a convolution integral of a jet function, J(x, kf , fcj), with the matrix element involv- 



4 Since the messenger modes from Ref. 18] do not spoil factorization for cases with a product of color 
singlet soft and collinear operators, we can neglect them in our analysis. 



ing the collinear fields, (tt\O c (x)\0) which gives (j)^(x), and that involving the soft fields, 
(E c (?/, s')\O s (kj~)\Ab(v, s)). In what follows we only need the form of the soft operator 

0.(k+) = [(h$S)tPL(SthV>)] [{dS) kt iP L {S^u) kt \ , (12) 

where h$' is an HQET heavy quark field, S is a soft Wilson line, and the subscripts denote 
the momentum carried by the fields. For our purposes the most important aspect of the 
analysis is that O s only involves the the combination h^^Pih^. Thus, by heavy quark 
symmetry 



[E c (v',s')\0 8 \A b (v,s)) = -=uv c (y',8') (Y ~ v ,fJ -)^P L u Ab (v,s) X,,, 



1 

(Z* c {v',s')\O s \A b (v,s)) = u^{v>,s')iP L u Ab {v,s)X^ (13) 



where v and v' are the four-velocities of the A& and £M respectively. The spinor field 
normalizations are u(v, s) u(v, s) = 1 for the A& and £ c , and u a (v, s) u a (v, s) = —1 for the 
£*. is the most general vector compatible with the symmetries of QCD, 

Xfj, = an M + bvp + cv'^. (14) 

Note that in Eq. ()12|) the part of O s involving the light quark fields is parity violating, 
so we need not worry about the fact that Aft — > S c is an "unnatural" transition. Using 
m Ab v = m-£ c v' + E n n to eliminate the term proportional to in Eq. (j!4j) . it is easy to see 
that any term in X^ proportional to v' does not contribute, so only remains. Hence 
the ratio of the rates for Aft — > S c 7r and A b — > S*7r are determined model independently at 
leading order in Aqcb/itiq and a s (mg), similar to the B° — > D^n case. We find 

l[ A A b - ^ = 2 + 0[A QCD /m Q , a s (m Q )] . (15) 
1 (Aft —> L c 7r) l j 

To evaluate the square of the matrix element in Eq. (|13[). we used the spin sums from 
Ref. jlj] for the various spin states. The explicit calculation shows that the rate to X* 
with \s'\ = 3/2 vanishes, as required by angular momentum conservation. 

A practical complication in testing this prediction is that the states decay to A c tt, 
and so both decay channels A b — > E(*)°7r° — > A c 7r~7r and A;, -> EW+tt" -> A r 7r°7r contain 
a 7T° that makes the reconstruction hard at hadron colliders. This can be circumvented by 
studying Aft — > E(*)°p° decays. In this case the final states are A c 7T~n + n~ . Decays to a vector 
meson are potentially more complicated due to the fact that "long-distance" contributions 
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can induce transverse polarizations at the same order in Aqcd/-E , tt- However, at leading 
order in a s (mg) these long-distance contributions vanish for the p° final state and we 
obtain 



T(A b - SOpO) * n 
It is also worth noting that 

Y(B° -> _D°7T ) r(A b -> S^vr ; 



(9 



AqcdMq , a s (m Q ) • (16) 



O 



r (£o ^ D o p o) r ( A6 ^ s o p o) 1 - L A Q CD/m « ' a ° ( v m Q A Q CD j J ' 
where in contrast to Eqs. (fT5|l and (fTfljl this prediction requires a perturbative expansion at 



the intermediate scale yAq CD mQ. 

The decays A& — > H C A" decays are also Cabibbo-allowed. (These decays involve u ss 
popping" so only E®K° is allowed, not 'E^K~). They are similar to A b — > S^7r in the sense 
that the leading contribution in the heavy quark limit vanishes. As shown in Table H] there 
are three S c "ground states", S c , H' c , and S*. The S c and E£ can mix, but the former is 
expected to be mostly the state that transforms as 3 under flavor SU(3), while the latter 
is mostly a 6. The H* also transforms as a 6, and forms a heavy quark spin symmetry 
doublet with the S' c . Thus, a relation similar to Eq. (|15j) also holds in this case, i.e., 
T(A fe -> H*iT)/r(A 6 -> H^) = 2 + C[A QCD /m Q , a s (m Q )]. This prediction may be hard 
to test since E' c decays to S c 7. One can also consider Cabibbo-suppressed A& decays, e.g., 
Af, — > H c 7r, and the weak decays of other baryons containing a heavy bottom quark. 

Perhaps the most exciting possibility is the existence of heavy baryonic pentaquark states. 
Recently several experiments claimed to observe a baryon + (154O) with the quantum 
numbers of K + n. A possible explanation is to consider the G + as a bound state of two spin- 
zero ud diquarks in a P-wave with an s antiquark If diquarks play an important role 
in these exotic states then the analogous heavy flavor states, C = c[u<i] 2 and G& = b [ud] 2 , 
may be below threshold for strong decays by AE ~ — 100 Me V and AE ~ — 160 Me V 
respectively fl9^. 5 Since the spin of the light degrees of freedom is si — 1, we expect from 
heavy quark symmetry that Gq come with a doublet partner of similar mass, Qq, as shown 
in Table HJ with a mass splitting of order Aq CD /mQ. From the mass splittings for the S c 
and S c we expect m®* — rriQ c ~ 70 MeV and me* — me b ~ 22 MeV. In this case the Qq may 



5 It is possible that the Qq arc above the strong decay thresholds ji^] . The assumptions in our analysis are 
that (i) Qq decay weakly; and (ii) the spin of the light degrees of freedom is si = 1, as suggested by [lsj . 
If (i) is correct but (ii) is not. it would be easy to modify our predictions, including Eq. l(T§|) . 
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FIG. 3: Comparison of the weak nonleptonic decays of the A& and 0;,. 



also be stable with respect to the strong interactions and decay to 0q7- Since the splitting 
for is larger, it is possible that the G* is just above the strong decay threshold, making 
the spectroscopy even more interesting (like in D* decays). 
It may be possible to discover the Qb,c via the decay chains 

0^e°7T + and 6° -> e+vr" -> K s pir~ -> ix + n-pix- (18) 

that are Cabibbo-allowed and lead to all charged final states. The most interesting aspect of 
the 0^ — > 0° decay is that in the diquark picture the correlation is maintained, as shown in 
Fig. |3J and so no additional suppression factor is expected. In weak 0^ decays to ordinary 
baryons this would not be the case. While we do not know the 0q production rates, we can 
estimate the branching ratios in Eq. (|18|). The lifetime of a weakly decaying ©b jC is expected 
to be comparable with other weakly decaying hadrons that contain a charm or a bottom 
quark. The 0^ — > 0°7r + amplitude factorizes, and is related to 0^~ — > Q^iv via a formula 
identical to Eq. (JHJ). For the nonleptonic rate we obtain in the heavy quark limit, 

r(ef->egir+) _ m| t (l-r|)» 1 W 4 r_ fu7 e ^(1 + 18r 2 + r 4 ) (19) 

F(A b - A c7 r") " <(1 - rlf C « ax ) 2 144r| I'N^J ^ + ^ + **> (19) 



4r/i(w® ax )r 72 « ax ) r (l - r|) 2 (l + r, 



2 (l-r 2 ) 4 



where vq = mQ c /mQ b , and 77^2 are the two Isgur-Wise functions that parameterize the weak 
©6 — > matrix elements where ^i(l) = 1. In particular 

(6 c (t/, s') I cr& |e 6 (v, s)) = ^ [^77! H - v a v'%(w)] u(v', s') 75(7* + v' a )T(j p + v p ) l5 u{v, s) . 

(20) 

Thus, B(Qfr —> ©c7r + ) is expected to be similar to B(A& — > A c n). If the Qq states exist 
then an analysis of the Aq CD /?71q corrections would be warranted, as the mass of the light 
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degrees of freedom is sizable. We expect £>(0° — > © + 7r~) to be at the few percent level, 
while the other branching ratios in Eq. (JTSj) may be of order unity. 

In summary, we studied nonleptonic A& decays to A c ir, S c 7r and S*7r. Eqs. (jlUj) . (|15|). 
and (|19|1 are our main results. In the mq ^> Aqcd limit the — > A c n rate is related 
to Afc — > A c £is, and we found that r(A& — > A c tt) is expected to be larger than T(B — > D^tt), 
as observed by CDF. At leading order in Aqcd/^q the Af, — > S^vr rates vanish, but an 
analysis of the leading contributions suppressed by Aq CD /mQ was still possible. We predict 
r(A fe -> S» /T(A 6 -> S c tt) = T(A b -> E*p) /r(A 6 ^ S c p) = 2 + 0[A QCD /m Q , a s {m Q )\. 
We also discussed properties of pentaquarks with a b or c, including a possible discovery 
channel if they decay weakly. 
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